Draft version February 5, 2008 

Preprint typeset using I^T^X style cmulatcapj v. 6/22/04 



O 
O 

(N 

Oh- 
< 

> 

(N 



O 

o 

Of 

6 

CO 

c3 



X 



COLLAPSE OF NEUTRON STARS TO BLACK HOLES IN BINARY SYSTEMS: A MODEL FOR SHORT 
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ABSTRACT 

The accretion of « 0.1 — 1 Mq of material by a neutron star through Roche lobe overflow of 
its companion or through white-dwarf/neutron-star coalescence in a low mass binary system could 
be enough to exceed the critical mass of a neutron star and trigger its collapse to a black hole, 
leading to the production of a short gamma-ray burst (SGRB). In this model, SGRBs would often 
be found in early-type galaxies or in globular cluster environments, though they could also be formed 
through stellar-wind accretion in high-mass binary systems of star-forming galaxies. Although the 
collapse event is likely to be of short (<C 1 s) duration, afterglow complexity could be produced by 
neutron decoupling in the fireball and subsequent accretion of the companion remnant, for example, 
the disrupted white-dwarf core. 

Subject headings: gamma-ray bursts — X-rays: binaries — stars: winds, outflows 



1. INTRODUCTION 

Recent discoveries from Swift and HETE-II have led 
to the identification of host galaxies for short GRBs 
(SGRBs), which is the class of GRBs with durations < 2 
s that tend to be spectral ly harder than the classic al long- 
duration GRBs (LGRB; iKouveliotou et aLlllQQfl. GRB 
05050 9B, the first rapidly localized SGRB l)Gehrels et all 
2005), displayed a single short spike with duration of 
40±4 ms. No radio or optical counterpart emission 
was detected. A large elliptical cD galaxy at redshift 
z = 0.225 overlapped the 9.3" error r adius of the XRT o n 
Swift, and is the likely host galaxy (Blo om et al.l l2006'l. 
No supernova emissions were detected to faint limits, rul- 
ing out an ori gin associated with supernovae, includin g 
Type la SNe ijHiorth et afll2005l iK^Ikarni et aU l2005l. 
Its apparent isotropic energy was w 4.5 x 10 48 ergs, with 
a smaller absolute energy release if the emission was col- 
limated. 

In contrast, GRB 050709, a SGRB discovered with 
HETE-II, is associated with a n irreg ular star-forming 
galaxy at z = 0.160 l)Fox et aU^OO^ . This burst dis- 
played a 70 ms single-spiked 30 - 400 keV hard pulse 
follow ed by a 130 s long smooth bump (jVillasenor et al.l 
|2005|) . The energy released in the bump was about the 
same as the energy released in the spike. X-ray flares 
occurred w 100 s and w 16 days after the GRB. The af- 
terglow was also detected at optical frequencies, with a 
jet break apparent at 10 days after the burst. A beam- 
ing factor /f, w .03 implies an a bsolute energy release of 
w 2 x 10 48 ergs (IFox et alJI200l . 

GRB 050724 is a Swift-detected SGRB that exhib- 
ited a 250 ms spike followed by an w 3 s 15 - 150 
keV hard bump and an « 100 s 15 - 25 keV soft bump 
starting 1 s and » 4 s after the trigger, respectively 
ijBarthelmv et alJ l2~005). This burst is clearly associated 
with an elliptical galaxy at z = 0.258. The X-ray af- 
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terglow of GRB 050724 is remarkable in that it exhibits 
a rapid decline in the X-ray flux at rj 200 - 400 s af- 
ter the trigger — similar to the behavior seen in LGRBs 
(Taghafer ri et al.ll2005|) — and an X-ray flare peaking w 
1/2 day later. A break in the optical light curve » 1 da y 
after the GRB suggests f b w 0.01 (jBerger et al J 1200^1 . 
implying an absolute total energy of « 4 x 10 48 ergs. 
The energy released in the soft bump and the later X- 
ray flare is comparable to the energy measured in the 
prompt phase. 

Besides these bursts, two other SGRBs have also been 
rapidly localized. GRB 050813 is associated with a 
cluster elliptical galaxy at either z = 0.722 or z > 1 
(JBergerl 120051 IGladders et all 120051. GRB 050906 is 
possibly associated wi th a nearby galaxy at z = 0.03 
(jLevan fc Tanvirl 12005 L The association of 3 out of 5 
SGRBs with elliptical galaxies excludes models associ- 
ated wit h massive stars as a co mmon origin fo r these 
sources (jProchaska et all 120061) . iGal-Yam etaD (|2005f) 
perform an archival analysis of SGRBS localized through 
the Interplanetary Network to enlarge the sample to in- 
clude two sources at z < 0.15 and two others at z > 0.3. 

Here we consider whether the collapse of a neutron star 
(NS) to a black hole (BH) i n a binary system, as con- 
sidered earlier for LGRBs bvlVietri fc Stella! l)1999j) and 
for SGRBs by iMacFadven et alJ l)2005|) . either through 
accretion-induced collapse (AIC) or through the coales- 
cence of the low-mass companion with the NS, can pro- 
duce SGRBs. Statistical arguments for the model are 
presented in Section 2, explanations for X-ray flares are 
considered in Section 3, and a discussion and summary 
are given in Section 4. 

2. STATISTICS OF NEUTRON STARS IN BINARIES 

Observations of low mass X-ray binaries (LMXBs) and 
high-mass XBs (HMXBs) in the Galaxy with both NS 
and BH companions provide support for the idea that 
the collapse of NSs to BHs could occur in binary sys- 
tems. The most widely considered evolutionary pathway 
to form a BH binary (BHB) involves instead the collapse 
of the evolved core of a massive (« 25 - 100 Mq) star in a 
primordial binary directly to a BH with mass « 10 M 
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(for review, see iLewin et al.lll995|) . Depending on the 
range of initial stellar masses that collapse to NSs and 
BHs, the rate of production of L MNSBs in the Galaxy 
is 10 ~ 5 A_R yr -1 , with A_5 ~ 1 (jPorteeies Zwart et al.l 
1997), at least by a factor «s 30 greater than the forma- 
tion rate of LMBHBs. This should lead to a preponder- 
ance of NSXBs over BHXBs in our Galaxy, contrary to 
observation. If the mean lifetime of NSXBs is smaller 
than BHXBs, then this discrepancy could be resolved. 
One possibility is that a large fraction of NSXBs in low- 
mass systems form BHXBs through AIC. In this case 
» 10 4 - 10 5 such events should occur over the lifetime 
of the Galaxy. Another possibility is that the number of 
NSXBs could be reduced relative to BHXBs if the rate 
of coalescence events, limited by the measured rate of 
SGRBs, is comparable to the LMB formation rate. 

Approximately 2 5% of BATSE GRBs are SGRBs 
l)Kouveliotou et all 119931) ■ implying that SGRBs occur 
at the rate 
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throughout t he universe fo r a full-sky BATSE GRB rate 
of 550 yr -1 {Band 2002), where /& is the mean beam- 
ing factor of SGRBs. Swift is less sensitive to SGRBs 
than BA TSE, and detects onl y 1 SGRB every 2-3 
months (|Barthelmv et al] 12005ft . or < 10% of the Swift 
GRBs. The present data are inadequate to determine 
the mean redshift z$grb of Swift SGRBs, which may, 
moreover, be different from the mean redshift of BATSE 
SGRBs. The value {V/V max ) SGRB = 0.39 ± 0.02 mea- 
suring the deviation from a Euclidean distribution for 
BATSE SGRBs is significantly larger than for BATSE 
LGR Bs ((V/V max ) LGRB = 0.29 ± 0.01 : rGuetta fc Piranl 
120051 and references therein), but is also sufficiently dif- 
ferent from 0.5 to suggest that cosmological effects are 
starting to play a role. 

The total number of LMXBs in galaxies of different 
types is poorly known, because LMXBs may be quies- 
cent for most of their lifetime. Measurements (j^ilfanov 
l2004j) show that, within a factor of w 2, the number of 
luminous (> 3 x 10 35 ergs s _1 ) LMXBs is proportional 
to the stellar mass of the host galaxy for spirals and el- 
lipticals. If the total number of LMXBs is proportional 
to the number of luminous LMXBs, then we can esti- 
mate the production rate Almx b of LMXBs in the local 
universe using the LMNS binary production rate in our 
G alaxy per unit stellar mass im plied by the calculation 
of iPorteeies Zwart et"aT] l|1997|) . The result for z < 1, 
using a stellar mass in the Milky Way of 3 x 1O 11 M0, is 
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where (« 0.1; see lRudnick et ai1 l2003'l is the fraction 
of baryon mass in the form of stars in the local universe, 
the ratio of baryon density compared to the critical den- 
sity is 5%, and we use a Hubble constant of 70 km s _1 
Mpc -1 . This rate exceeds the SGRB rate, eq. (Q, if 
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where e (< 1) is the fraction of NSs in LMBs that collapse 
to a BHs over the ~ 10 10 yr lifetime of galaxies. Unless 
A_5 3> 1, a significant fraction of NSs in binaries must 



therefore collapse to BHs to provide the observed rate of 
SGRBs. 

A contact or short-period LMXB forms through spiral- 
in of the lower initial mass secondary during a common 
envelope phase while the higher initial mass primary is 
evolving through a giant phase. If the primary evolves 
to form a NS before the secondary coalesces with the he- 
lium core of the primary (which depends on the initial 
separation and initial masses of the stars), then a close 
LMNS binary will form. When the secondary evolves 
through its hydrog en-burning phase and expa nds to over- 
fill its Roche lobe l)Podsi adlowsk i"e"t al.ll20 02'l , additional 
accretion onto the NS can trigger collapse to a BH and 
the formation of a SGRB. Alternately, the secondary can 
lose angular momentum through gravitational radiation 
and magnetic braking to spiral in and coalesce with the 
NS. If the pulsar had evaporated the companion, all that 
may remain of the secondary is a helium core. 

Formation of a BH through AIC or white-dwarf/NS 
coalescence requires that the NS accretes enough mass 
to exceed the maximum stable mass of a NS. De- 
pending on the equation of state of nuclear matter , 
this mass is w 2.2 - 2.9M (|Kalogera fc Bavmi ri996;i. 
though rotation could incre ase these values by < 20% 
l|Shapiro fc Teukolskvlll983^ . Sufficient mass could also 
be accreted from stellar winds to trigger AIC in HMBs. 
These systems would be found only in star-forming galax- 
ies, and may display distinctive signatures due to the 
different environment. 

Recent observations show that NSs with masses signif- 
icantly exceeding the Chandrasekhar mass are found in 
LMBs. PSR J0751+1807 is a 2.1 M© millisecond pulsar 
in a 0.26 d o rbit around an » 0.2M© helium white dwarf 
companion ijNice et alJl2005|) . As its orbit decays, the 
white dwarf will overflow its Roche lobe. A number of 
evolutionary outcomes are possible, including an ultra- 
compact X-r ay binary phase and tidal disruption of the 
white dwarf ijErgma et al.lll997j) . Given the large mass 
of the NS, the NS may also collapse to a BH when it ex- 
ceeds its maximum mass through accretion. NS masses 
> 1.7Mq are found for seve ral millisecond puls ars in 
the globular cluster Terzan 2 ijRansom et alJ l2005'). The 
masses of most millisecond pulsars are, however, near the 
Chandresakhar limit, even allowing for th e mass needed 
to spin up the pulsar ijJacobv et alJl2005|) . but selection 
biases may exist against the detection of high mass NSs 
in millisecond pulsar systems if magnetic field decay ac- 
companies accretion. 

The total mass accreted by the compact companion 
in LMXBs over the lifetime of the Galaxy is p s 10 5 M^, 
using the LMXB energy output computed by iGil fanov 
(|2004t eq. (19)), for a 10% accretion radiation efficiency 
and a bolometric-correction factor of « 2. Most of this 
energy is radiated by sources with luminosity <2x 10 38 
ergs s , the Eddington limit of a 1.4 M object. If NSs 
in ~ 10 5 LMNSBs each accrete about the same fraction 
of this mass during that time, then a large fraction will 
collapse to BHs in the Galaxy. The proportionality of 
the number and luminosities of LMXBs in systems with 
old stellar populations compared to the Milky Way in- 
dicates that this mechanism will also operate effectively 
in early-type galaxies. The larger fraction of LMXBs in 
globular clusters, which are more prevalent in elliptical 
galaxies ijSarazin et alJ l20031 , and the long time required 
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for AIC, means that SGRBs would often be found at 
large distances from the center of E/SO galaxies. 

3. X-RAY FLARES FROM SGRBS 

One of the surprises of the Swift and HETE-II results 
is the complex X-ray afterglow behavior of SGRBs. It is 
hard to invoke density inhomogeneities for bright flares 
observed hundreds of seconds to days or more after the 
trigger, especially in those cases where SGRBs are asso- 
ciated with old stellar populations that should inhabit a 
rarefied quasi-uniform ISM. 

There are, however, at least four other processes that 
can produce episodes of emission from SGRBs on dif- 
ferent timescales. First is the duration of engine activ- 
ity. Before reaching the moment of collapse, the NS in 
a LMXB could have accreted enough mass and angular 
momentum to have been spun up to millisecond peri- 
ods and would, in effect, be a recycled pulsar. As in 
the case of a hyp ermassive NS formed in binary neutron- 
star coalescence ijShibata et al.l l20061. the newly formed 
black hole would be surrounded by a dense torus that ac- 
cretes on subsecond timescales, accounting for the first 
burst of emission of the SGRB. Millisecond structure in 
the lightcurve could be explained by angular-momentum 
transport through the magneto-rotational instability in 
the material accreting on the newly formed BH. 

A second timescale that could appear in the X-ray 
lightcurve of SGRBs relates to the deceleration of the 
blast wave as it sweeps up material from the ISM. For 
an apparent isotropic energy release of Eq — 10 50 £ , 5o 
ergs from a blast wave with mean entropy per baryon 
To = 300r3oo that decelerates in a uniform surrounding 
medium with density n = 10 _2 n_2 cm -3 , the decelera- 
tion timescale 
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The value of td represents the time when unabsorbed syn- 
chrotron emission from a decelerating blast wave reaches 
its maximum value. The value of Tq is unknown for 
SGRBs, so the afterglow synchrotron emission could ex- 
plain emission peaks on various timescales, 3 for example, 
emission bumps observed a few seconds to minutes after 
the start of GRBs 050709 and 050724. The peaking of 
this emission could however be concealed by the decay 
phase of the prompt emission if td is very short, which 
could happen if Tq 3> 300. 

Another emission episode can occur because of neutron 
decoupling from the fireball during its expansion phase 
Iperishev et alJ Il99flt iBahcall fc Meszarosl 12000ft . The 
baryon load would be rich in neutrons from the accret- 
ing neutron-star material. If the elastic nuclear interac- 
tion timescale is longer than the comoving time, neutron 
decoupling occurs when To/Tdcp = u (> 1), where the 
decoupling Lorentz factor 
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3 If SGRBs are sources of 100 MeV - GeV radiation, then 
GLAST will be able to constrain lower limits of To from 77 trans- 
parency arguments. 



Here a pn = 30 mb, y is the ratio of neutrons to protons 
in the expanding fireball, and A/c is the characteristic 
timescale during which energy is injected, so that A/c = 
10 9 A 9 /c = 33A 9 ms. 

The Lorentz factor of neutrons when they decouple 
is r„ = Tow -4 / 3 , and the fraction of total energy in 
the decoupled neutrons is (1 + u 4 ' 3 /y) _1 . The neu- 
trons will decay before colliding with the forward proton 
shell when the comoving time t' n = x co i/T n c 3> 900 s, 
the mean neutron lifetime. The collision radius where 
the shell of neutron-decay protons impact the deceler- 
ating forward proton shell is x co i = 2 1 / 3 u 8 / 9 Xd, where 
x d = T 2 ct d /(1 + z) = 2.6 x 10 16 ( J B 5 o/n- 2 r| 00 ) 1 / 3 
the deceleration radius. If Xd ~3> 6 x 10 15 r 30 ou 
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cm, 



then most neutrons will decay before the impact. 

From this, we derive the observer time of the collision, 
given by t co i = 2 1//3 u 32 / 9 £ [ j. When u = 2, a neutron- 
decay proton shell with significant amount of energy will 
collide with the forward shell, resulting in a strong shock 
that could be observed as the second peak near 100 s 
in GRB 050724. After the shock crosses the shell and 
injection ceases, the flux rapidly declines. This results 
from fast synchrotron cooling and expansion of the shell. 
Note that only one flare can be produced from this pro- 
cess, thoug h other neutron de coupling scenarios can be 
considered l)Beloborodovll2003[) . 

The preceding three mechanisms may also equally op- 
erate for NS/NS and NS/BH mergers. In the case of 
AIC of a NS to a BH, however, the impact of gravita- 
tional and neutrino radiation produced at the time of 
collapse of the neutron star to a BH would heat and dis- 
rupt the companion star, causing the ejection of material 
outside the Roche surface. Fallback of this material onto 
newly formed black hole can trigger flares hours or days 
later, as sugg ested by simulatio ns of white dwarf merg- 
ers with BHs llFrver et al.lll999)) . The stellar companion 
can also be disrupted through tidal forces if its orbit is 
perturbed. Only » 1O _5 M0 of matter is required for 
X-ray flares from SGRBs, and fallback of the disrupted 
stellar material could make flares weeks to months after 
the event. We expect that a large fraction of BHs in 
LMBHBs are just above the maximum stable mass of a 
n eutron star if this model is correct. 

iMacFadven et all l)2005|) also argue that shock interac- 
tion of the ejecta with a red giant companion could make 
flares seen minutes after the prompt emission. Such a 
process is unlikely to explain the rapid decline in GRB 
050724 due to light travel-time delays from different parts 
of the red-giant surface. Moreover, because comparable 
X-ray energy is detected in the prompt emission and in 
the delayed X-ray flares, considerably more energy must 
be directed toward the red giant than is emitted in the 
beamed prompt emission. X-ray flares unassociated with 
SGRBs would be therefore be predicted from the un- 
beamed shock-heated red giant surface. 

4. DISCUSSION AND SUMMARY 

Although we have considered a model of SGRBs from 
the collapse of NSs to BHs in binarie s, as first pro- 
posed for SGRBs bv IMacFadven et all Ipfll . SGRBs 
may represent a heterogeneous class. The currently 
favored scenario for SGRBs is compact-object coa- 
lescence involving NS/NS and NS/BH binaries (brief 
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revie ws are given bv [Zhang fe Meszarosl 120041: iPirarJ 
120051 see iRuffert & Janka| J2fl01|) for Newtonian calcu- 
lations and iShibataetldT pOOa) for fully general rela- 
tivistic calculations). Such systems that will coalesce 
within a Hubble time are known to exist, and can re- 
lease sufficient ene rgy through vv — ^ e + e~ processes 
l) Janka et al.lll999t iLee et al l 12005]) to power a SGRB. 
Compact object merger scenarios with short delay times 
may, however, conflict with the statistics of SGRB host 
galaxy types and the properties of obser ved NS binaries 
l)Gal-Yam et alJl2005t iNakar et al J 12006(1 . Evolutionary 
calculations of the timescale for AIC of NSs in LMXBs 
are required to determine if the model proposed here is 
in accord with the preference of SGRBs for early-type 
galaxies. 

The features that make compact-object coalescence 
an attractive model for SGRBs also apply here, be- 
cause both involve bursting events derived from the grav- 
itational energy released by the collapse of a supra- 
massive NS to a BH. Initial loss of angular momen- 
tum support along the rotation axis of the NS would 
form evacuated axial regions and a massive, dense 
torus of neut ron-star material that would accrete o n 
ms timescales i|Vietri fc Stellall998tlShibata et al.l2006jl . 
Along the poles, a baryon- dilute jet would be cre- 
ated through vv annihilatio n <|Alov et alJl2005H or MHD 
Iprenkhahn fc Spruit! [2002) processes. Insofar as the 
progenitor systems would not involve stellar winds or su- 
pernovae, the collimated outflows would be launched into 
a low-density external medium. 

The occurrence of flaring emissions long after the 
SGRB takes place was not, however, predicted in NS/NS 
or NS/BH coalescence scenarios. For NS/NS merg- 
ers, the collapse i s over after only a few tens of ms 
l) Janka et aU 11999). The torus formed by a transient 



highly magnetized hypermassive NS collapses to a BH 
formed in a NS/N S merger lasts for only ~ 10 ms 
(ShibaTa~eFaIl2006). Even in the case of NS/BH coales- 
cence where the NS is shredded, acti vity is not expecte d 
to last for more than a few seconds l)Davies et al.l l2005). 

In the model proposed here, the disruption of the com- 
panion star during coalescence could produce flares long 
after the SGRB. In the case of AIC, by contrast, if the 
stellar companion could survive the collapse of the NS 
to a BH and the accompanying impact of neutrinos and 
gravitational radiation, then the system may relax to a 
BHB system. Depending on the type of companion star 
at the time of collapse, either a LMBHB or HMBHB will 
be produced. This scenario offers an alternate pathway 
for BHB and BH microquasar formation. 

Gravitational radiation is emitted in this model 
through two mechanisms. One is the collapse of the NS 
to the BH. The collapse of deformed, rapidly rotating 
NSs with millisecond periods would enhance the gravita- 
tional wave output compared to collapse events that are 
nearly spherically symmetric. The coalescence of white 
dwarf cores with masses of a few tenths of a Mq onto NSs 
could also be detectable, and we suggest that template 
waveforms for such events be calculated for comparison 
with LIGO data. 
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